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ABSTRACT 

The presence of superconducting and superfluid components in the core of ma¬ 
ture neutron stars calls for the rethinking of a number of key magnetohydrodynamical 
notions like resistivity, the induction equation, magnetic energy and flux-freezing. Us¬ 
ing a multi-fluid magnetohydrodynamics formalism, we investigate how the magnetic 
field evolution is modified when neutron star matter is composed of superfluid neu¬ 
trons, type-II superconducting protons and relativistic electrons. As an application of 
this framework, we derive an induction equation where the resistive coupling origi¬ 
nates from the mutual friction between the electrons and the vortex/fluxtube arrays 
of the neutron and proton condensates. The resulting induction equation allows the 
identification of two timescales that are significantly different from those of standard 
magnetohydrodynamics. The astrophysical implications of these results are briefly dis¬ 
cussed. 
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1 INTRODUCTION 

Existing at the extremes of physics, neutron stars serve as 
excellent cosmic laboratories. Some of the most striking fea¬ 
tures are related to their magnetic properties. Measured 
magnetic field strengths, generally inferred from the star’s 
dipole spin-down, by far exceed the strengths of terrestrial 
magnets. Observations also suggest a link between the vari¬ 
ous classes of neutron stars. Old millisecond pulsars thought 
to be formed in low-mass X-ray binaries have fields between 
10 s — 10 10 G, while the fields of ‘classical’ rotation-powered 
pulsars range between 10 10 — 10 13 G. A third class of slow- 
rotating, highly magnetised neutron stars, so-called mag- 
netars, reaches field strengths up to 10 15 G. This class is 
believed to include both, soft gamma repeaters and anoma¬ 
lous X-ray pulsars. Understanding the long-term evolution 
of the stars’ magnetic fields might be key to establishing con¬ 
nections between the different cla sses and forming a unified 
picture of the neutr on star ‘zoo ’ (iKaspii 120106: IVigano et alj 
l2013l:lHardin3l2013h . 

Unsurprisingly, such enormous strengths suggest that 
magnetic fields are crucial for the neutron stars’ d ynam¬ 
ics. A s first pointed out by Thompson & Duncan (1 19951: 

1 19961) . magnetic field decay on a timescale of ~ 10 4 yr could 

* E-mail: vanessa.graber@soton.ac.uk 


power the high activity of magnetars. The rotational en¬ 
ergy is not sufficient to explain the observed emission of 
these objects. There are also observations indicating that 
the magnetic dipole fields of s tandard pulsars evolve on a 
timescale of the order 10 7 yr (Lvne. Manchester fc Tavlod 
1 19851 : iNaravan fe Ostriker] 1 199C ). The actual mechanisms, 
causing the magnetic field to change on these rather 
short timescales, are only poorly understood and there 
is no definitive answer to the question of which part of 
the neutron star dominates the magnetic field evolution. 
Most theoretical studies and numerical simulations focus 
on the crust as the s ource of the field decay and neglect 
the core contribution j|Pons_& 1 Ge]3 perj|2007l : IVigano et al.l 
l2013l : lc!ourgouliatos fe Cummindl2014l) . However, one could 
argue that the core, which carries the majority of the star’s 
inertia and magnetic energy, should also play a role in the 
magnetic field evolution. 

The problem of magnetic field evolution in isolated 
neutron stars has been discuss ed by a number of authors. 
iGoldreich fe Reiseneggeil ll 19921 ) determined several mecha¬ 
nisms that are present in an ionised plasma consisting of 
neutrons, protons and electrons. Ohmic diffusion due to 
the interaction of relativistic electrons and lattice nuclei 
causes magnetic field dissipation in the crust. This mech¬ 
anism is most effective on small scales and, thus, not ex¬ 
pected to affect the large scale evolution of the crustal field. 
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However, Ohmic decay could be enhanced by Hall drift, 
which is in itself conservative but may redistribute mag¬ 
netic energy from large to gradually smaller lengthscales. 
The combined effect, sometimes referred to as the Hall cas¬ 
cade, could cause field evolution on a timescale of order 10' 
years JGoldreic.h fc Reiseneggei]ll99 Z). However, recent neu¬ 
tron star crust simulations show no strong cascading be¬ 
haviour but suggest the existence of a quasi-equilibrium es- 
tablished on timescales shorter than the Ohmic timescale 
llPons fe Gepperj 2010l: Gourgouliatos fe Cummind 12014 


IWood fe Hollerbach 20151 ). In the core, standard Ohmic 
decay is negligible because the interior is expected to 
form a type-II superconductor (iBavm. Pethick fe Pines! 
Il969bl . see also below). Electron-proton scattering, al- 
ready acting on very long time scales in the star’s interior 
(IBavm. Pethick fe Pines! Il969al ). is then restricted to the 
normal conducting cores of fluxtubes. These only contribute 
a tiny fraction to the star’s total cross section. Hence, the 
coupling timescale is increased further, making this dissi¬ 
pation mechanism irrelevant. Additionally, ambipolar diffu¬ 
sion, describing the motion of charged particles and mag¬ 
netic field lines relative to the neutrons, could cause mag¬ 
netic field decay and drive the flux from the core to the 
crust (Goldrcich & Rciscncg gerlll992l). This mechani sm was 
originally considered by Thompson fe Duncanl ll 19951 ) to ex¬ 
plain the magnetar activity. However, more recent results 
seem to indicate that the timescale for ambipolar diffusion 
considerably increases when the superfluid nature of the 
neutrons or proton superconductivity is t aken into account 
(iGlampedakis, Jones fe Samuelssonll201ll ). 

Lacking a clear answer, the question of magnetic field 
evolution in a neutron star is revisited in this paper. We 
foc us on the outer core and use the formalism developed 
by_ Glam^edakiSj^ndersson^Samuelsson (201l]) (see also 
iMendell fe Lindblom 1991 . Mendelll ll99ll ). which presents 
a general set of macroscopic hydrodynamic equations for 
a multi-fluid mixture. We address, in particular, the effect 
the superconducting component has on the evolution of the 
magnetic field. 

The presence of superfluid and superconducting com¬ 
ponents in neutron stars is firmly supported by observa¬ 
tions and microphysical calculations. Traditionally, glitches 
and post-glitch relaxation timescales on the order of months 
to years are seen as observ ational evidence of superfluid¬ 
ity. lAnderson fc Itohl (11975! ) first proposed that the neu¬ 
tron star’s dynamical evolution during and after a glitch 
could be explained by the weak viscous properties of a su¬ 
perfluid component that is coupled to the crust. Moreover, 
recent spectral analyses of the neutron star in the super¬ 
nova remnant Cassiopeia A indicate that the surface tem¬ 
perature of this young object decreases faste r than one 
woul d expect from s tanda rd cooling models (IPage et al] 
l201ll ; IShternin et al.l 1201 ll ). The rapid cooling could be 
explained by enhanced neutrino emission, resulting from 
the onset of neutron superfluidity and proton supercon¬ 
ductivity in the core. In addition to observations, the¬ 
oretical calculations provide strong reasons for macro¬ 
scopic quantum condensates in neutron star cor es. The 
i dea o f superfluid interiors was first put forward by iMigdaj 
T959I). several years before the first detection of a pulsar 


( Hewish et al .tjH). A few hundred years after birth, neu¬ 


tron stars are in thermal equilibrium and have tempera¬ 


tures of 10 s — 10 s K |Tsuruta|l998^ Page, Geppert fe Webeil 


120061 : iHo, Glampedakis fe Anderssonl 201 2I ). While this is 
certainly hot in terms of terrestrial physics, the temper¬ 
atures lie well below the Fermi te mperature o f nuclear 
matter, which is of the order 10 12 K ( Sauldfl989l ). Apply- 
ing the microscopic theory of labo ratory superconductors 
jBardeen, Cooper fe SchriefTerll 19571) to the neutron star in¬ 
terior suggests that the neutrons and protons form Cooper 
pairs and, thus, condense into a superfluid and a supercon¬ 
ductor, respectively. 

The properties of the proton superconduc- 
tor were further discussed in the seminal paper by 
IBavm. Pethick fe Pinesl (Il969bl ). As the conductivity of nor¬ 
mal co nducting matter is very large dBavm. Pethick fe Pinesl 
Il969al ). the authors argued that flux could not be expelled 
from the star’s interior. The phase transition into a super¬ 
conducting state has to occur at constant magnetic flux. 
Two characteristic lengthscales determine the state of the 
superconducting protons. The London penetration depth 
describes the magnetic field’s exponential fall-off from the 
surface of a superconductor or the fluxtube cores due to the 
Meissner effect. The coherence length, which characterises 
the typical dimensions of a Cooper pair, is equivalent to the 
core radius of a fluxtube. The ratio of these two quantities, 
the Ginzburg-Landau parameter, k, dictates the type of su¬ 
perconductivity. For k > l/\/2, it is energetically favourable 
for the material to increase the surface area between normal 
and superconducting regions, implying that magnetic 


flux can penetrate the medium in the form of q u antised 
fluxlines (fAbrikosovlll957l) . IBavm, Pethick fe Pinesl (1 1969bl ) 
estimated k for the neutron star’s outer core and found 
that the protons would enter a metastable type-II state by 
forming a regular array of fluxtubes (see als o Figure 1 of 
IGlampedakis. Andersson fe Samuelssonll201ll ). 

The presence of fluxtubes significantly influences the 
magnetic properties of the star because the flux is no longer 
locked to the charged plasma but is mainly confined inside 
the fluxtubes. As standard coupling mechanisms, like Ohmic 
dissipation, are suppressed as a result of pairing, interac¬ 
tions of fluxlines with their surroundings determine the mag¬ 
netic field evolution on macroscopic lengthscales. The most 
prominent of these effective coupling processes, known as 
mutual friction, is the scattering of electrons off the fluxtube 
magnetic hel d dAlpar, Langc r fe Saulsl Il984 iMendelll Il99ll : 

I Andersson. Siderv fc Qomerll200f)l ). In the following, we de¬ 
rive an equation for the magnetic field evolution of a su¬ 
perfluid/superconducting mixture using a smooth-averaged 
formalism. We translate the mesoscopic phenomena, influ¬ 
encing individual fluxtubes, to the large scale picture by 
applying the framework known from standard resistive mag¬ 
netohydrodynamics (MHD). 

The structure of this paper is as follows. In Section [2] 
we introduce the hydrodynamical equations for a multi-fluid 
system and the corresponding Maxwell equations. In Section 
[3] the standard resistive MHD formalism is reviewed. We 
then discuss the case of superconducting matter in Section 
H focussing on the standard r esistive coupling introduced by 
Alpar, Langer fe Saulsl (11984 ) . After giving the most general 
form of the superconducting induction equation, we use a 
few assumptions to simplify and interpret our results. We 
finally conclude with a discussion in Section [5] 

Note that we will work in an inertial frame and carry 
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out the analysis in a coordinate basis. Vectors are, thus, de¬ 
noted by their individual components. We also use Gaussian 
units in the remainder of this paper. 


2 THEORETICAL BACKGROUND 
2.1 Multi-fluid hydrodynamics 

The model present ed in this section largely builds upon 
the re cent work by iGlampedakis, Andersson fc Samuelssonl 
l|201ll) . Making use of the Lagrang ian formalism developed 
by Carter, Prix and collaborators ([Carter fe Langloiall995l : 
IPrbd 120041 : lAndersson fe Comer! l2006h . a full set of MHD 
equations for the superfluid/superconducting bulk in the 
outer neutron star core is derived. The simplest represen¬ 
tation of this Fermi liquid is a mixture of three components, 
namely relativistic electrons, superconducting protons and 
superfluid neutrons. In the following, the constituents are 
denoted by roman indices, x = {e,p,n}. Note that in 
order to keep the discussion clear, we neglect the pres¬ 
ence of muons. However, generalising to the four-constituent 
case would be straightforward as electrons and muons are 
strongly coupled a nd move as on e component on macro¬ 
scopic lengthscales JMendelllll99ll '). 

The dynamics of the three-fluid system is governed 
by two Euler equations, one representing the superfluid neu¬ 
trons and the other representing the electron-proton con¬ 
glomerate. The charged fluids can be characterised as a 
single component if charge neutrality holds over macro¬ 
scopic distances, i.e. n e = n p , where n x denotes the par¬ 
ticle number density. As the electrons are mobile enough 
to quickly equilibrate any local charge imbalances, this re¬ 
quirement is fulfilled in the neutron star core (j.Tacksonl 
Il999l : iGlampedakis. Andersson fc Samuelssonl 12011 1. Addi¬ 
tionally, the electron mass is significantly smaller than the 
proton mass, m e m p , which allows us to neglect any elec¬ 
tron inertial terms. The resulting macroscopic Euler equa¬ 
tions are 

\dt + UnVjj (t)‘ + EnW pn ) - E n Wp n V’r“ 

= —V* (pn + *3?) + /mf + ,/mag,n 3 (1) 


+ UpVj] (v p + £pWnp) + E p wi p X7 l V? 

= -V i (A + 3>)--/m f +/ia g ,p- (2) 

P P 


The averaged fluid velocities are denoted by v £ and the mass 
densities by p x = rnn x , where m = m n = m p is the baryon 
mass. wj y = uj — Uy is the relative velocity, >1? the gravi¬ 
tational potential and £ x the entrainment parameters. By 
definition, the latter satisfy the condition n p e p = n n e n . The 
specific chemical potentials are defined as 


~ _ /^n 

Mn = 5 

m 


, — Mp + Pe 
[l — 


(3) 


The two Euler equations are supplemented by three conti¬ 
nuity equations for the number densities, 

d t n x + Vi (n x v^J = 0, (4) 


species, and the Poisson equation 

V 2 <f> = 4 t rGp. (5) 


p = p x is the total mass density of the system and G the 
gravitational constant. 

The variational approach used to derive the Euler 
and continuity equations explicitly distinguishes between 
the fluid momenta and velocities. This formalism provides 
the possibility to include any changes, caused by the su¬ 
perfluid and superconducting condensates, into the hydro- 
dynamical model. In contras t to the mome ntum equations 
of standard plasma physics (|jacksonlll999l l. Equations JlJ 
and m incorporate new inertial terms due to entrainment, 
which arise from the strong coupling of Fermi liquids, and 
terms that go beyond the standard electromagnetic interac¬ 
tion given by the Lorentz force. For the fluid mixture in the 
outer neutron star core, the right-hand sides of the momen¬ 
tum equations contain the total magnetic and mutual fric¬ 
tion forces per unit volume, /mag.x and /mf, respectively. The 
former one is caused by interactions of the vortex/fluxtube 
magnetic field with the charged fluid. We point out that the 
neutron fluid experiences this magnetic force because pro¬ 
tons are entrained around each neutron vortex and create 
an effective magnetic field. In the absence of entrainment, 
the magnetic force on the neutron component would vanish. 
Finally, the mutual friction forces arise from the dissipa¬ 
tive coupling of the vortex and fluxtube array with the fluid 
components. 

We keep in mind that our hydrodynamic model is 
solely based on averaged quantities and reflects the macro¬ 
scopic behaviour of the fluid components. It is on these large 
scales that we have a method to deal with the presence of 
the quantum condensates in a consistent way. Taking advan¬ 
tage of the large numbers of vortices/fluxtubes, we average 
over the respective arrays and obtain a smooth-averaged de¬ 
scription of the magnetic and mutual friction forces. Using 
this formalism, it is possible to determine how the presence 
of vortices and fluxtubes influences the macroscopic dynam¬ 
ics of a neutron star, i.e. its rotational and magnetic evo¬ 
lution. If individual vortices/fluxtubes do not overlap and 
are distant enough so that interactions within one array can 
be neglected, the averaging procedure is obtained from the 
macroscopic qu anti sation conditions originally developed by 
iOnsagerl Il949l l and lFevnmanl (Il955l l for the dynamics of ro¬ 
tating superfluid helium. Assuming that neutron vortices 
and proton fluxtubes are locally straight and directed along 
the unit vectors, itj and k p , the arrays can be assigned vor¬ 
tex and fluxtube surface densities, A/” n and N p , respectively. 
As the vorticities, W x , are related to the circulation of the 
averaged canonical momenta, the macroscopic quantisation 
conditions are given by 

vv‘ = Vj ( Vk + £n wl n ) = NnE n , (6) 

Wp = (vl + e p wl p ) + a p B l = N p k p , (7) 


where 4 = kk x points along the local vortex direction with 
the quantum of circulation 

k = « 2.0 x 10 -3 cm 2 s“\ (8) 

2m 

and we define 


cip — 


e 


9.6 x 10 3 C“ 1 s" 1 . 


which reflect the conservation of mass for each individual 


me 


(9) 
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The proton charge, the speed of light and the Planck con¬ 
stant are denoted by e, c and h, respectively. 


2.2 Macroscopic magnetic induction 


In the averaged framework, the total magnetic induction, 
B l , is the sum of three individual components, namely the 
averaged fluxtube and vortex field and the London field, 

B i = B l p + Bi + bl. (10) 

The former two contributions are obtained by multiplying 
the surface densities, A/” x , with the flux carried by a sin¬ 
gle line of the lattice, </> x . The fluxes can be derived by 
considering the dynamics of a single vortex/fluxtube on 
mesoscopic lengthscales, denoted by bars on the respective 
quantities. Using the corresponding quantisation condition 
and the mesoscopic Ampere law (e.g. see Ap pendix Al of 
iGlampedakis. Andersson fo Samuelssonl (120111 ) for details), 
it is possible to derive generalised London equations for the 
mesoscopic magnetic fields, B j), 

A*V 2 S* - B x = -(p x K^8(f), (11) 


where 5(f ) is the two-dimensional delta function located at 
the centre of each vortex/fluxtube, </ x is defined below and 
the effective London penetration depth is given by 


A» 


f 1 1 £n £ p 

\ ^Pp a p 1 — £n 


1/2 


( 12 ) 


In the absence of entrainment, e n = e p = 0, this expres- 
sion reduces to the standard result of superconductivity 
jTinkhaml 120041 '). Taking advantage of the symmetry and 
using cylindrical coordinates, the inhomogeneous Helmholtz 
equation m can be solv ed using a Green’s functi on ap¬ 
proach in two dimensions dFetter fc Hohenberdll969l h Inte¬ 
grating the resulting magnetic induction over a disc of radius 
r>A, perpendicular to K* gives for the magnetic flux 

J Bids = <£ x k’- (13) 


For a proton fluxtube, we obtain the expected unit of flux, 


<f>p = 4>o = — = — ~ 2.1 x 10 ' Gem 2 , 
a p 2e 

whereas the flux of a superfluid vortex is 


(14) 


1 £n 


4 > 0 - 


(15) 


We note that the minus sign originates from it), and B l n 
pointing into opposite directions. In the absence of entrain¬ 
ment, the neutron flux would be zero. The averaged contri¬ 
butions from the two arrays to the macroscopic magnetic 
induction, B', are then obtained by 

Bi = A4</> x /t x - (16) 


The third contribution to the induction is the London 
field. It is a fundamental pro perty of a supercondu ctor and 
associated with its rotation (iTillev fe Tillevl Il99tf ). While 
superfluids need to form vortices in order to support any 
circulation, the superconducting fluxtube array is not re¬ 
lated to the macroscopic rotation. However, these dynamics 
induce an additional magnetic field inside the superconduc¬ 
tor, whose axis is parallel to the rotation axis. In contrast to 


the London field is not of microscopic origin but related 
to the macroscopic electromagnetic current (see Subsection 
im Combining the quantisation conditions © and 0 with 
Equation m, the London field can be related to the macro¬ 
scopic fluid properties. Assuming that the hydrodynamical 
lengthscales are sufficiently small to ensure constant entrain¬ 
ment parameters, we have 

bl = 1 ~ £n ~ £p e ijk V jV l (IT) 

dp -L S n 

The proton entrainment parameter is related to the effective 
proton mass via 

ml , , 

Sp — 1-£ « 0.3, (18) 

m 

where we have used the es timate m* « (0.6 —0.9)m given by 
IChamel fc Haensel (|2006l ~) for the outer neutron star core. 
For small proton fractions, £ p = p P /p <S 1, an approx¬ 
imation which is valid in the interior of a neutron star, 
the neutron entrainment coefficient is negligible because 
e n ~ XpEp -C 1. Taking the proton fluid to be tightly coupled 
to the neutron star’s crust through the magnetic field and, 
thus, rotating rigidly at the observable pulsar frequency, we 
can substitute a canonical rotation period to calculate an 
estimate for the magnitude of the London field. Using the 
normalisation Pio = P/( 10 ms), we find 

b L « 9.2 x 10 _2 .P/q 1 G, (19) 

which is many orders of magnitude smaller than the mag¬ 
netic field strengths usually invoked for neutron star physics. 
Hence, it is generally justified to neglect the London field in 
Equation Goj, an approach we will take in Subsection ED 
We can simplify the magnetic induction, B x , one step 
further by taking the properties of the vortex and flux¬ 
tube arrays into account. Although the individual fluxes, i/> x , 
are comparable, the contribution from the superconducting 
protons dominates: Consider smooth-averaged fluid veloci¬ 
ties of the form e yfc V jVk = 2 Q x z l . In this case, the vor- 
tices/fluxtubes are aligned with the z-axis, i.e. z l = £*, and 
Equation © gives 

•A/n = — [Hn + e n H pn ] ■ (20) 

As we expect the lag, H pn = , to be small on macro¬ 

scopic scales and the neutrons to be coupled to the crust, 
the vortex surface density is 

Af n ~ « 6.3 x 10 s P^ 1 cm" 2 , (21) 

The neutron vortex density is, thus, fixed by the rotation of 
the neutron star. Using previous estimates for the entrain¬ 
ment parameters, the magnetic field strength of the neutron 
vortex array is 

B n « JV n EpcPo « 4.0 x 10“ 2 P 1 ( ) 1 G. (22) 

This is again many orders of magnitude smaller than typi¬ 
cal neutron star field strengths, which implies that the mag¬ 
netic field of the outer neutron star core is mainly confined 
to the proton fluxtube cores. It is, therefore, important to 
investigate which mechanisms affect the motion of individ¬ 
ual fluxtubes in order to link the small scale behaviour to 
the large scale evolution of the star’s magnetic field. The 
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fluxtube density can be estimated to 

jV p = ^«|-»i4.8x 10 18 B \2 cm" 2 , (23) 

00 00 

with the normalised magnetic field B 12 = B/( 10 12 G). 


2.3 Macroscopic Maxwell equations 


In order to capture the electromagnetic response of the fluid 
mixture correctly, the Euler equations m and 0 have to 
be supplemented by Maxwell’s equations. Taking these to 
be valid in our multi-fluid mixture, we have to redefine, or 
rather reinterpret, the various fields accordingly in order to 
make Maxwell’s equations suitable for a type-II supercon¬ 
ductor. 


As mentioned before, the London held, despite be¬ 
ing of small magnitude, plays an important role for the 
electro dyna mics. As discussed by ICarter, Prix fe Lan gloisl 
d2000h and iGlampedakis. Andersson fe Samuelsson ( 201 lh . 
the London held is closely connected to the macroscopic 
electromagnetic current, 


T 


CTleWp e , 


(24) 


which enters the macroscopic Ampere law. In contrast to 
standard MHD, where the equality H l = B l is satished, 
the averaged magnetic induction, B l , and the macroscopic 
magnetic held, /7\ are no longer equivalent in a type-II su¬ 
perconducting sample. Instead, Ampere’s law reads 

e ijk VjH k = e ijk Vjbl = — J\ (25) 

c 

where the displacement current has been neglected because 
the fluid motion is slow compared to the speed of light. This 
deviation from standard MHD can be also understood in 
terms of the classihcation generally app lie d to terrestrial 
superconductors (see for example lTinkhamll2004h . In labo¬ 
ratory experiments, one distinguishes between macroscopic 
electromagnetic currents that generate a macroscopic held, 
H 1 , and magnetisation currents only affecting the meso¬ 
scopic induction, which is B l in our notation. A supercur¬ 
rent, circulating around each vortex/huxtube and generat¬ 
ing is attributed to the second class. It does not con¬ 
tribute to the held TP = 6^, which is created by the current 
J 1 . Hence, the macroscopic magnetic induction, B l , given in 
Equation m differs from the magnetic held, Ft 1 . For com¬ 
parison, in vacuum or normal conductors, no magnetisation 
currents are present and the identiheation H' — B l = B l 
can be made. In the present case, 

e ijk X7jH k = e ijk VjB k =* ^ j. (26) 

In addition to Ampere’s law, we use 


ViB i = 0, 


(27) 


which has to hold everywhere in the superconducting Huid, 
and the macroscopic Faraday law, 

dtB 1 = -ct ijk S7jE k . (28) 


Instead of dehning the macroscopic electric held as the av¬ 
erage over the microscopic equivalent, we take advantage of 
the remaining huid degree of freedom, namely the electron 


Euler equation, to obtain an expression for E l . Neglecting 
again the electron inertial terms, we have 

, 1 iik er , m e „i, ^ FI , , 

c e ca p p p 

where FI represents the total force exerted on the electrons 
due to interactions with the surrounding huid components. 

Combining Equations (1281) and (1291) leads to an evolu¬ 
tion equation for the magnetic induction that only depends 
on macroscopic fluid variables. However, the procedure relies 
on the forces, FI, and we address this in the following sec¬ 
tions. 


3 MAGNETIC FIELD EVOLUTION IN 
STANDARD MHD 

Before discussing the more complicated problem of magnetic 
held evolution in a superhuid/superconducting mixture, we 
briefly review the approach taken in normal resistive mat¬ 
ter, which allows us to compare our new results with a well 
studied model. 


3.1 The MHD induction equation 


In a charged plasma containing electrons and protons, the 
only relative flow present is the motion between the two 
components. Assuming that a frictional mechanism would 
damp these dynamics and try to bring the two species into 
co-motion is straightforward. Hence, resistive coupling act¬ 
ing on a timescale r e leads to a dissipative force on the elec¬ 
tron fluid, 


F i = 

e 


n e m e 



(30) 


Substituting this force into Equation (1291) gives a generalised 
Ohm’s law 


E i = -- e ijk 


J en e 


m T 1 

5/C-—V l (£e + $) + — (31) 

e (Je 


CPpCLpCTe 

— 


with the standard electrical conductivity defined by 

(32) 
leads to an 


-X7 l B m 


Equation m combined with Faraday’s law 
evolution equation for the magnetic induction, 


dtB 1 = 


(' v l pB m ) 


ijk-r-, 

-e J V j CfcZn 


47r<J e 


iik 

e 


V j €.klm 


WIC s p 
47rep p 


C V S B P ) B m 


(33) 


We have used Ampere’s law for normal matter (1261) to elim¬ 
inate the macroscopic current in the last expression. The 
second and the third term on the right-hand side represent 
the Ohmic decay and the Hall evolution. We can extract the 
following well-known timescales, 


and 


47ra e L 2 

TOhm = q 
C 2 


4nep p L 2 

mcB 


(34) 


( 35 ) 
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where L is the characteristic lengthscale over which the mag¬ 
netic field changes. 


3.2 Flux-freezing in MHD 


We can estimate the two charact eristic timescales for a 
neutron star core. According to iBavm. Pethick fe Pines! 
(Il969al l. the electrical conductivity associated with the inter¬ 
action of highly relativistic electrons and normal, degenerate 
protons is given by 

a e « 5.5 x 10 28 T s - 2 p% 2 (^ 5) 3/2 s-\ (36) 

where Tg = T/(10 8 K) is the star’s normalised tempera¬ 
ture, P 14 = p/( 10 14 g cm -3 ) the normalised total density 
and a; p the proton fraction. Approximating the characteris¬ 
tic lengthscale by the radius of the neutron star, the Ohmic 
diffusion timescale is 

rohm « 2.4 x 10 13 T 8 ~ 2 Ll p% 2 yr, (37) 

with the normalised lengthscale Lg = L/(10 6 cm). For the 
Hall timescale, we obtain 

THaii W 1.9 x 10 10 B ±2 Ll Pl4 (^) yr. (38) 

Both estimates are many orders of magnitude larger than the 
typical spin-down ages of radio pulsars. We would, therefore, 
expect Ohmic decay and Hall term to be negligible for the 
evolution of the plasma’s magnetic field. In this idealised 
case, which is commonly used to approximate astrophysical 
or laboratory plasmas, the induction equation reduces to 

dtB i = e yfc V i e Wm («p£ m ) • (39) 

Using Equation imp, we can rewrite the last expression and 
simplify the result using the standard Lie derivative, 

d t B i + C Vp B i = - B'Vjvi. (40) 


The left-hand side describes how the magnetic field vector, 
B l , is transported with the fluid flow, v p . Taking into ac¬ 
count that the mass of the proton plasma is conserved, we 
use Equation to further simplify, 


d_ 

dt 



+ -Cu. 



= 0 . 


(41) 


This implies that the magnetic field is moving with the fluid, 
i.e. the fluxlines are frozen into the proton plasma. 

As soon as Ohmic and Hall terms play a role for 
the dynamics, this frozen-in condition is destroyed and field 
lines are no longer forced to follow the protons. In particu¬ 
lar, if Ohmic decay characterised by the conductivity, cr e , is 
included, the induction equation resembles a diffusion equa¬ 
tion. It encodes how the magnetic field lines diffuse through 
the fluid and reconnect, leading to the decay of magnetic en¬ 
ergy as discussed in Section [3.31 If the Hall term is present 
but Ohmic decay is negligible, the induction equation re¬ 
duces to 

d t B i = e ijk V j ekim [v l e B m ^j . (42) 

In contrast to Equation (1391) , the electron velocity enters the 
magnetic evolution law. This implies that the relative mo¬ 
tion between electrons and protons becomes important and 


the magnetic field is frozen into the electron fluid. The Hall 
term in Equation (1331) itself is not dissipative but may act 
to redistribute magnetic energy from large scales to smaller 
ones, where it can decay ohmically. Many studies of the in¬ 
duction equation’s non-linear b ehaviour are based on results 
from hydrodynamic turbulence (lKolmoeorovlll94~th , as it has 
sever al similarities with the vorticity equation of a viscous 
fluid (IGoldreich fe Reiseneggeii 119921 '). However, recent nu¬ 
merical simulations in the context of neutron stars have 
shown no evidence of strong cascading behaviour. Instead 
the Ha ll cascade appears to be saturated at long length- 
scale s (iPons fe Gepper tl 120101; iGourgouliatos fe Cummind 


120141 ; [Wood fe Hollerbach 2015lh 


3.3 Magnetic energy in standard MHD 


The conservative and dissipative nature of the different 
pieces in Equation (1331) is illustrated by considering the evo¬ 
lution of the magnetic energy. In order to compare the stan¬ 
dard MHD plasma with the superconducting mixture later 
on, we calculate the magnetic energy associated with the 
work done by the Lorentz force. In its standard form, the 
Lorentz force density, given by 




B j V j B i - {B k B k ^j 


(43) 


is composed of a tension and a pressure term. The work is 
obtained by calculating the product with the position vector, 
r 1 , and integrating over the volume, V. Using the product 
rule and Equation m , we arrive at 


W L = J nFt dV=^-J V* (rjB^Bi - ^B 2 r?j dV 


1 

47T 


B'B.VV, - -B 2 X7 i n dV. 


(44) 


The total gradient term can be rewritten as a surface integral 
using Gauss’ theorem. As no discontinuities are present at 
the boundary of the plasma region, we can push the integra¬ 
tion radius to infinity. Provided that the magnetic induction 
vanishes at infinity, the surface contribution is zero. The in¬ 
tegrand of the second term in Equation m simplifies to 
the well-known magnetic energy density 

52 dU = J £ mag dV. 


W L = 


87 r 


(45) 


Changes in the magnetic energy are, thus, determined by 


d£ u 


dt 


d_ 

dt 


Bj dB 1 
47t dt 


(46) 


Calculating the product of the induction equation with Bi 
and using the product rule to rewrite the result, we find 

2 


d£ rnag _ 1 Jsptf 7 D 
dt _ Hr £ (Vs " pj 


ik vlB K - 


4na t 


t ijk {S7 J B k 


me 

ditepp 


e ijk e jlrn (\' l B rn )B k 


V‘E ;. 


(47) 


The last term contains all contributions that can be writ¬ 
ten as a divergence. After integrating over the volume, it 
is possible to convert this part into a surface integral using 
Gauss’ theorem. Additionally, the third term has to be zero 
due to the properties of the Levi-Civita tensor. As expected, 
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the Hall term is conservative and does not contribute to 
the change in the magnetic energy density. Using Ampere’s 
law (1261) and the generalised Ohm’s law (13111 , the remaining 
terms are simplified to 

dSmas 1 T* ? 7->fc 't 

». — -J tijk VpB - 

8t C P < 7 e 

- V* [dESi - — (Ae + $) Ji 1 , (48) 

L47T e J 

where S z = N k EjBk is the Poynting vector. Equation ([481) 
clearly shows that any resistive plasma is subject to the 
decay of magnetic energy due to Ohmic diffusion and the 
energy loss is proportional to J 2 . The inertial term vanishes 
if the protons are not able to move, which is, for example, 
the case in a standard metal, or the macroscopic current and 
the proton velocity are aligned. 


4 MAGNETIC FIELD EVOLUTION IN 

SUPERCONDUCTING NEUTRON STARS 

4.1 The coupling force: ‘standard’ resistivity 

We now return to the question of magnetic field evolution 
in the superconducting outer neutron star core. In order 
to apply a formalism similar to the resistive MHD dis¬ 
cussion, we need to determine the forces, FI, exerted on 
the electron component by the various fluid constituents 
and the vortices/fluxtubes. However, this is where things 
get complicated. Due to the multi-fluid nature of the 
superfluid/superconducting mixture, there are not simply 
two components coupled by a single resistive force. We 
could imagine a variety of ways for the components to 
inter act with eac h other ranging from electron scatter- 


1 ° 1—1 

19841: Andersson. Siderv 

& Comer 2006() and vortex- 

fluxtube interactions 

Ruderman. Zhu & Chenl 19981: 

Jahan-Miri 

2000|; Linkl 

20031) to shear or bulk vis- 

cosity ( 

Andersson. Comer & Glampedakis 20051: 

IShternin & Yakovlev 2008: iManuel. Tarrus & Tolos: 2013if. 


Choosing a more pedagogical approach to our problem, 
we pick one specific mechanism, determine how it affects 
the electrons on mesoscopic scales and translate this into 
a macroscopic picture. While we won’t provide a complete 
picture of the magnetic field evolution in the core, this 
method provides more insight to how different mechanisms 
could play a role. 

We keep in mind that the magnetic field is locked to 
the superconducting fluxtubes and their motion determines 
the evolution of the magnetic field. We, therefore, consider 
the scattering of electrons off the vortex/fluxtube magnetic 
fields as a source of mutual friction. This ‘standard’ resis¬ 
tive coupling in a superfluid/superconducti ng mixture, first 
discussed by lAlpar, Langer fe Sauls! dl984ll , results in two 
forces acting on the electrons, 


F 1 = F l -I- F l 

e pe i ne 


: F x 

pe* 


(49) 


We neglect the contribution from electrons scattering off 
the neutron vortices because N p A/” n . This implies that 
electron-fluxtube interactions are markedly more common 
and, thus, dominate the electron coupling. The macroscopic 
force, Fl, is obtained by multiplying the electron drag force, 


fl, exerted on a single fluxtube, with the fluxtube density, 

N p , 

Fp e =N P fl= UpppFTL (vl - Up) . (50) 


The dimensionless drag coefficient, 1Z, contains all the infor¬ 
mation about the coupling on mesoscopic scales and u p is 
the velocity of a single fluxtube. 

In order to determine an evolution equation for the 
macroscopic magnetic field in superconducting matter, we 
have to eliminate any quantities from Equation (1501) that are 
defined on mesoscopic lengthscales. Hence, the next step is 
to rewrite the fluxtube velocity, u p , in terms of the macro¬ 
scopic fluid variables. This can be achieved by using the force 
ba lance for an individu al fluxtube, an approach introduced 
bv [Hall fe VinerJ (|l956i ) for the description of superfluid he¬ 
lium. We have 


F — fit + /m + ft + fl m — 0, (51) 


where the fluxline inertia is neglected. Our force balance 
equation includes the electron drag force given above, the 
Magnus force, /^, the tension force, fl, and the electromag¬ 
netic Loren tz force, fl m . The d i fferent forces have b e en cal - 
culated by iGlampedakis, Andersson fe Samuelssonl (l201ll l 
and are given by 


fl = 


H c \K 
4-7T a D 


e ijk k p {vl-ul), 

(52) 

k J pVjk l p, 

(53) 


where H c i is the low er critical field for superconductivity 
dTillev fe Tillevl[i990l 'l. and 


/ ■% ink * p p 

em = Pp^e 


(54) 


Note at this point that we are interested in the linear anal¬ 
ysis of one specific resistive mechanism. For this reason, our 
force balance EH does not include a ‘pinning’ force, result¬ 
ing from t he magnetic short-range interaction between the 
two arrays jRuderman, Zhu fe Chen||l998|:| Jahan-MirHl2000l : 
lLinkll2003l : IGlampedakis fe Anderssonll201lib 

Calculating repeated cross products of the force bal¬ 
ance equation with k l p , pointing along the local orientation 
of a fluxtube, it is possible to express the mesoscopic flux¬ 
tube velocity in terms of the averaged fluid velocities, 


Up = vl + - + 1 ^ 2 (iZfl + e ljk k p f^ , (55) 

where 

/: = e ijk k»wl p + — (fi+ fL) ■ (56) 

PpK \ / 


Combining the previous relations, we observe that the first 
term in Equation (1561) and the electromagnetic force, /* m , 
cancel each other. Then, the effective force, fl, is equivalent 
to the fluxtube tension, 


fl = - ft = 

PpK 


Ft cl „ n _, ^ i 

- - HitWjKp. 

47ra p p p 


(57) 


Substituting Equations (ESJ) and E3 back into Equation 
ESI finally gives for the macroscopic drag force, 

H c \(j>oNptZ 


Ft = - 


4tt(1 +1Z 2 ) 


(tZk J p VjKp + N k k p k l p Vik p ^j . (58) 


For a straight fluxtube array, the tension force and, thus, 
the electron coupling would vanish. 
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4.2 The superconducting induction equation 


Having determined the force, FI, exerted on the electron 
component due to scattering off the fluxtube magnetic fields, 
we use Equation for the macroscopic electric field to 
derive a generalised Ohm’s law that is valid in the super¬ 
fluid/superconducting mixture, 


E i 


1 ijk 

— t 

C 


VjBk - (/ie + 4>) + 


-Hcl^pA/'p 

47rca P M P 


X I +K2 (^pVjKp + . (59) 

With Faraday’s law (12811 . we obtain a superconducting in¬ 
duction equation describing the evolution of the macroscopic 
magnetic field in the outer neutron star core, 


d t B i 


e ijk Vj 


€-klm 



Hci 00-A/p 

47ra P M P 


tz 


l+TZ 2 


(UKpVlkl + £fcIm«p«pV s Kpj 


(60) 


Let us specify the l ower critical fi eld, H c i, to simplify this 
expression further dTillev fc Tillevfl99f)l '). The field is related 
to the energy per unit length of the fluxtube, £ p , via 


H cl 


An£ p 

<t> o 


(61) 


The fluxtube energy, on the other hand, is determined by 
the characteristic lengthscales of the superconducting phase 
and given by 


£ P — 



(62) 


Here, £ p denotes the proton coherence length and the ef¬ 
fective London penetr ation dept h was defined in Equa¬ 
tion (1121) . According to iTinkhaml (l2004lb one can estimate 
In (A*/£ p ) sw 2, which leads to 


p p n 


1 £n 


ppK m 
2n rrip 


27T 1 £n £p 

Combining the previous equations, we find 

K(j> oA/” p m 


dtB i = 


n 


CfcZr 


(v l e B m ) 


2n mi 


l+TZ 2 


(riKpVlkl + tfclmApApV s Kpj 


(63) 


(64) 


which is the main result of our paper. 


4.3 A simplified set of equations and the field 
evolution timescales 

At this point, it seems natural to make several assumptions 
about the actual physics of the multi-fluid mixture inside a 
neutron star in order to find a simplified version of Equation 
As discussed previously, the main contribution to the 
macroscopic magnetic induction is given by the fluxtubes. 
In this case, we can neglect the weak London field and the 
superconducting Ampere law (1251) dictates that the protons 
and electrons are co-moving on large scales, i.e. u p ~ vl, 
and the macroscopic current vanishes. This also implies that 


the local direction of the fluxtube array is aligned with the 
direction of the magnetic field because 

B i = B& « jV p (f) 0 Kp gives £*««;’. (65) 

Using these simplifications to rewrite the force on the elec¬ 
tron fluid, we obtain 

K ~ + e ijk BjB l V t B k ) . (66) 

The induction equation, on the other hand, reduces to 


d t B i 


ijkr-r 

e Vj 



kB m 

27r m p 


n 


l+TZ 2 


[TZB l ViB k + t klm B l B s V s B r ‘ 


(67) 


We will compare this form of the superconducting induction 
equation with the standard MHD result in Equation (1331) . 
As in the resistive MHD case, we can extract two timescales, 


Tclif 


2nL 2 (l + H 2 ) m* 


kR, 


_p 

m 


and 


ri _ 2nL 2 (l + TZ 2 ) m* 
TZ kTZ 2 m ’ 


( 68 ) 


(69) 


where L is again the characteristic lengthscale over which 
the magnetic field changes. The naming convention of the 
two timescales might seem arbitrary at this point, but our 
choice will become clear later on. 

We can estimate these two timescales provided the 
strength of the mutual friction is known. A method to 
calculate the dimensionless drag parameter, TZ, for the 
coupling of rela tivistic electrons and a singl e fluxtubes 
is discussed in Alpar . Lang er fe Sauls dl984l) (see also 


Sauls. Stein fe Serenel Il982l : lAndersson. Siderv fe Cornell 


2006t1 1 l . The authors give the relaxation timescale, r pe , for 


the ‘resistive’ interaction and include effects caused by the fi¬ 
nite fluxtube size and the increase in moment of inertia due 
to the coupling of electrons and protons on much shorter 
timescales. The relaxation timescale is related to the drag 
coefficient via 


TZ = (KjVpTpe) -1 , (70) 

which leads to the following numerical estimate 



This gives R < 1 and implies that the standard friction 
mechanism is rather weak. Adopting this limit, we can ap¬ 
proximate for the neutron star core, 

r diss « 3.1 x 10 11 Le p( 4 1/6 (^g) ~ 1/6 yr (72) 

and 

rcons « 1.3 x 10 15 Ll Pl - 4 2/6 (^g) ’ 2/6 yr. (73) 

We have used Equation Gi to estimate the effective proton 
mass, i.e. the entrainment parameter. 


1 Note that according to l.Tone 3 te00fih . magnetic scattering off in¬ 
dividual fluxtubes is suppressed for very large fluxtube densities. 
Instead electron scattering by a cluster of fluxtubes dominates 
the coupling, leading to a much smaller drag coefficient, 7 Z. 
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4.4 Flux-freezing and magnetic energy 

For conventional electron-fluxtube coupling, the timescales 
for the magnetic field evolution are rather long and the dy¬ 
namics of the macroscopic induction are dominated by the 
inertial term in the induction equation. In the weak mu¬ 
tual friction limit, we are, thus, left with an equation that is 
equivalent to the ones discussed in Subsection l3.2l The mag¬ 
netic field in the superconducting sample is frozen to proton 
fluid, which implies that the superconducting fluxtubes are 
locked to the proton plasma, i.e. Vp ~ Up. Hence, electrons, 
protons and fluxtubes are comoving on large scales, which 
is different to the weakly resistive case of standard MHD, 
where the relative motion between the charged particles was 
important. 

In order to determine whether the additional terms 
in Equation 63 are conservative or dissipative and which 
timescale dominates, we discuss the evolution of the 
superconducting magnetic energy. As before, we evalu¬ 
ate the energy associated with the work done by the 
magnetic force. However, in a superfluid/superconducting 
mixture the standard Lorentz force (S3 has to be 
changed accordingly. For a non-rotating star in the ab- 
senc e of entrainment the t o tal magnetic force is given 
bv dEasson fe Pethickl Il977l: lAkgfin fe Wassermanl 120081: 


iGlampedakis. Andersson fc Samuelsson 201ll: Landerll201 


F l = — 
mag 4 tt 


BN :i Ki - V* PpB 


, dH c 

dpp 


(74) 


where H l cl = . This expression also contains a tension 

and a pressure term but both scale with Hh and B 1 instead 
of B 2 . The work associated with this force is given by 


ww = - 


rjH^Bi- ppB^-n 


dV 


- J (hUb^u - PpB^VV^ dV, (75) 

where we have used the product rule and Equation {27}. 
Similar to the standard MHD case, the first term can be 
rewritten using Gauss’ theorem. However, in the supercon¬ 
ducting outer core, this contribution does not simply vanish 
because discontinuities are likely to be present at the fluids’ 
boundaries. The dynamics that might arise due the presence 
of a current sheet at the crust-core interface or the type- 
II to type-I transition region in the neutrons star’s inner 
core are only poorly understood and significantly complicate 
the problem. Incorporating these interfaces would require a 
much more detailed understanding of the microphysics in¬ 
volved. In the following, we, therefore, omit a discussion of 
the surface terms and focus on the much simpler problem of 
magnetic field evolution in the bulk fluid. 

Taking into account that H c i is a function linear in 
Pp, the second integral in Equation (1751 reduces to 


MHD given in Equation (1461 shows that the superconduct¬ 
ing nature of the mixture gives rise to a new contribution for 
the change in energy density. In contrast to resistive MHD, 
the evolution of matter and the magnetic induction are no 
longer decoupled in the condensate. Equation 113 demon¬ 
strates that modifying the properties of the superconductor, 
such as the lower critical field, H c i, alters the magnetic en¬ 
ergy. This implies that an evolving matter configuration can 
be closely linked to a changing magnetic field. 

The second term in Equation {TT} is similar to the 
result for normal conducting matter. Calculating the prod¬ 
uct of the induction equation (EZ1) with Bi and using the 
product rule to simplify, we find 

B^=e isp (VsB P ) 


kvlB K - 


kBTZ 2 


27t(1 + 1Z 2 ) m 


m B l S7iBi 


kBTZ m 
2tv(1 + 1Z 2 ) m* 


e ijk B j B l ViB k 


V*Ej. 


(78) 


As before, E® denotes the divergence terms. This equation 
bears some resemblance with the result {47} found in stan¬ 
dard MHD and we would equivalently expect to obtain a 
conservative and a dissipative contribution. In order to de¬ 
termine which of the terms are nonzero or zero, we rewrite 
the tension using the following identity 

B l ViBi = t ijk t ilm {ViB m )B k . (79) 


The second term in Equation (1781) is, thus, proportional to 
e isp {X7 s B p )B l X7iBi = e i3p {X7 s B p ) e ijk e jlm {ViB m )B k . (80) 


Analogous to the Hall term of standard resistive MHD, this 
vanishes due to the properties of the antisymmetric Levi- 
Civita tensor. Thus, the second term in the superconducting 
induction equation {67} is conservative and does not modify 
the total magnetic energy of the superconducting mixture. 
The last term in Equation (1781) . on the other hand, is pro¬ 
portional to 

e isp (V s B p ) e ijk B j B l ViB k = J*e ijk B j e klm JiB m , (81) 


where we defined the vector 

T = e ijk \/jB k . (82) 


Rewriting the remaining two Levi-Civita tensors in terms of 
Kronecker deltas, 8], gives the projection 

J i e ijk B j e klrn JiB m = J t J i - ( J l B (j j B^j . (83) 

Decomposing the vector J l into a component parallel to Bi 
and one perpendicular to the magnetic field direction, i.e. 

= J\\Bi + Jl, we see that Equation (1551) only depends 
on the component of J l that is perpendicular to Bi, 

- (j'bI) ( J j B ,-) = Jl. (84) 


ITmag.bulk = J ^J- dV = f £ mag.scdF. (76) 

Taking the time derivative of the magnetic energy density 
gives two contributions 

<9£ m a g ,sc B dH cl , Hd 5 dB i 

dt 2 tt dt + 2 tt i dt ' 1 ’ 

Comparison with the corresponding expression of standard 


Similar to the Ohmic term in standard MHD, we also retain 
a dissipative contribution to the total magnetic energy in 
the case of superconducting MHD. It is given by 


d£u 


dt 


B_dH ci 

2-7T dt 


i ^ cl -r* ,. 
+ 


k v\ 


\B k 


H cl uBTZ m 2 
Atv 2 {1 + TZ 2 ) m* J± 2tt 4 


( 85 ) 
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Having calculated the change in the magnetic en¬ 
ergy density, we can associate the timescale raiss, given 
in Equation ([72)1, with a dissipative mechanism. Compar¬ 
ing the numerical estimate to the Ohmic decay timescale 
Gil), we observe that the resistive coupling in a super¬ 
fluid/superconducting mixture acts on a timescale, which 
is two orders of magnitude smaller than the standard MHD 
diffusion, 

T diss 1r . — 2r r '2 —5/3 / \ 5 ^ 3 foct\ 

- -« 1.3x10 T 8 Pi 4 (—) ■ (86) 

TOhm VU.U5/ 

On the other hand, the timescale (1381) for the Hall evolution 
in a normal conducting plasma can be compared to r cons in 
Equation (1731) . In contrast to the standard MHD case, the 
conservative timescale emerging from the superconducting 
induction equation is several orders of magnitude larger than 
the Hall timescale, 


^«6.8xl0 4 H 12 pI 4 4 / 3 (-^-)' 4/3 . (87) 

We also note that due to the dependence on the dimension¬ 
less drag coefficient, 1Z, the dissipative term in the supercon¬ 
ducting induction equation governs the magnetic field evolu¬ 
tion, whereas in standard MHD the conservative Hall term 
acts on shorter timescales. In the latter case, the order of 
the two timescales is necessary for any cascading behaviour 
to take place; the Hall term drives the magnetic field to 
shorter lengthscales, where it can decay ohmically. However, 
if diffusion is dominating the evolution, the redistribution of 
magnetic energy will happen on much longer timescale not 
causing a cascade. Hence, despite the close similarities be¬ 
tween the conservative terms in Equation m and standard 
M HD, we conjecture that the ana lysis of the Hall cascade 
bv lGoldreich fc Reiseneggerl h 1992h is not transferable to the 
superconducting case. 


5 DISCUSSION 

Strong magnetic fields are a key ingredient for many phe¬ 
nomena observed in neutron stars. Understanding the fields’ 
long-term evolution might give insight into the ‘metamor¬ 
phosis’ between the different neutron star classes, the chang¬ 
ing fields of standard radio pulsars or the high activity of 
magnetars. As the mechanisms causing field changes are only 
poorly understood, we revisited the question of magnetic 
field evolution and, in particular, discussed the influence of 
a superconducting component. Our aim was to rethink key 
notions of magnetohydrodynamics and develop a better in¬ 
tuition for the magnetic field evolution in a superconductor. 

In this work, we used a multi-fluid formalism to de¬ 
scribe the mixture in the outer neutron star core. The model 
introduced in Section [2] translates the presence of meso¬ 
scopic vortices/fluxtubes into the large scale dynamics of 
the fluid. As an application of this framework, we anal¬ 
ysed the conventional dissipative mechanism, i.e. the scat¬ 
tering of electrons off the fluxtube magnetic field. Based 
on the approach of standard resistive MHD, we combined 
a generalised Ohm’s law with Faraday’s law and the re¬ 
spective force on the electron fluid to derive a supercon¬ 
ducting induction equation. Considering the London field 
as a negligible contribution led to a simplified equation, 


which should be applicable to most astrophysical scenar¬ 
ios. Caution is in order when discussing highly magnetised 
objects. For field strengths above the upper critical field, 
B > H c 2 ~ 10 16 G dTillev fc Tillevlll990l ~). the superconduct¬ 
ing state breaks down and our averaged formalism no longer 
applies. According to iGoldreich fe Reiseneggeil (I1992T ). am- 
bipolar diffusion could potentially become important in this 
regime and drive field decay on the order of typical magnetar 
ages. 

To compare our new results for magnetic field evo¬ 
lution with the standard MHD case, the magnetic energies 
associated with the total magnetic forces were calculated. 
In our analysis, we significantly simplified the problem by 
omitting a detailed discussion of the surface terms, a key 
problem which needs to be addressed in future studies. This 
implies that effects originating at the crust-core interface or 
the type-II to type-I transition in the inner core, which could 
potentially drive the magnetic field evolution, are not taken 
into account. Instead, we focused on the evolution of the av¬ 
eraged magnetic field in the bulk. We found that in the limit 
of weak mutual friction, the inertial term dominates the field 
evolution. The fluxtubes move with the proton fluid and the 
flux is, as in the standard MHD case, frozen to the charged 
particles. We additionally showed that the new induction 
equation contains a dissipative and a conservative contribu¬ 
tion, similar to the Ohmic and the Hall term in normal con¬ 
ducting matter. However, the evolution timescales extracted 
from the superconducting induction equation for weak mu¬ 
tual friction, 10 11 yr and 10 15 yr respectively, are notably 
longer than the typical spin-down ages of neutron stars. 
We conclude that the conventional mutual friction mech¬ 
anism cannot serve as an explanation for the field changes 
in pulsars or the activity of magnetars, which would require 
timescales of order 10 7 yr and 10 4 yr, respectively. Simply in¬ 
creasing the strength of the mutual friction cannot provide 
a solution to this problem either. Due to the ^-dependence 
of Tdiss, the minimum dissipation timescale one could obtain 
with a frictional mechanism of the form m is 


_ 4-7tL 2 m p 

7m in — 

k m 


1.4 x 


io 8 l 2 6P - 4 1/6 



yr 


( 88 ) 


for 1Z = 1. We note at this point that all numerical estimates 
crucially depend on the lengthscale L. It can be identified 
with the curvature radius of the magnetic field and we chose 
the neutron star radius, R, to normalise the previous results. 
Recent work on field equili bria in superconducting neutron 
star cores by [Lander] (120141 ') suggests that the field configu¬ 
ration actually supports structures on a shorter lengthscale 
of L m 10 s cm. Adopting such an estimate would reduce the 
characteristic timescales by two orders of magnitude. In par¬ 
ticular, the minimum dissipation timescale, r m in, would be 
shortened to a million years, which is closer to the timescales 
of astrophysical interest. 

Our results additionally suggest that the highly con¬ 
ducting neutron star core might affect the crustal field and 
slow down its evolution. Making a precise statement at this 
point is, however, not possible due to the poorly known 
physics at the crust-core boundary. This transition is crucial 
in understanding how changes of the core magnetic field are 
translated to the crust. The analysis presented in this paper 
does, therefore, not reconc ile the discrepancy between short 
crustal decay timescales (IPons, Miralles fe Geppertl l2009l l 
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and the much longer core evolution. In order to significantly 
reduce the latter different dissipative mechanisms have to 
be invoked. The typical candidate for strong coupling is 
vortex-fluxtube ‘pinning’ due to the short-range magnetic 
interaction between the two arrays. While we did not ad¬ 
dress pinning specifically, discussing the vortex-fluxtube in¬ 
teraction would be the natural continuation of this paper 
as our prescription can deal with any coupling mechanism. 
Based on a mesoscopic description of the pinning process, 
one would have to determine how the coupling affects the 
electron fluid and substitute the respective force, FI, into 
the generalised Ohm’s law. Determining the superconduct¬ 
ing induction equation that would result from the pinning 
interaction will be left for future work. 
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